Quantum dynamics of a resonator driven by a superconducting single-electron 
transistor: a solid-state analogue of the micromaser 



D.A. Rodrigues, J. Imbers and A.D. Armour 

School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, United Kingdom 

We investigate the behavior of a quantum resonator coupled to a superconducting single-electron 
transistor tuned to the Josephson quasiparticle resonance and show that the dynamics is similar 
in many ways to that found in a micromaser. Coupling to the SSET can drive the resonator into 
non-classical states of self-sustained oscillation via either continuous or discontinuous transitions. 
Increasing the coupling further leads to a sequence of transitions and regions of multistability. 
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Systems where a mesoscopic conductor such as a 
single-electron transistor is coupled to a nanomechani- 
cal resonator have been studied intensively because the 
current through the conductor can be extremely sensitive 
to the motion of the resonator and hence may be used to 
monitor its position with almost quantum-limited preci- 
sion [l|, 0, H, 0] ■ Furthermore, where either the coupHng 
between the electrons and the resonator is non-linear [5| 
or the electronic transport occurs via a resonance [1] , dy- 
namic instabilities in the resonator can occur leading to 
self-sustained oscillations. The way a nanomechanical 
resonator can be driven into states of finite amplitude 
oscillation by successive interactions with a current of 
electrons in a conductor parallels the behavior of quan- 
tum optical systems, such as the micromaser, in which an 
electromagnetic cavity is pumped by interactions with a 
steady stream of individual two-level atoms Q . This con- 
trasts with a standard laser (a nanomechanical version of 
which was envisioned in where an oscillator interacts 
simultaneously with many two-level systems. 

In a superconducting single-electron transistor (SSET) 
transport can occur via resonant processes involving both 
coherent motion of Cooper pairs and incoherent quasi- 
particle tunneling, the simplest of which is the Joseph- 
son quasiparticle (JQP) resonance 8]. In the vicinity 
of a JQP resonance, the dynamics of a resonator cou- 
pled linearly to the SSET is very sensitive to the bias 
point 0, 0, • For bias points on one side of the res- 
onance, the SSET acts on the resonator like a thermal 
bath and its current can monitor the position of the res- 
onator with exquisite sensitivity. In contrast, biasing on 
the opposite side of the JQP resonance can drive the res- 
onator into states of self-sustained oscillation ^ . 

In this Letter we explore the quantum dynamics of a 
resonator coupled to a SSET and show that it is analo- 
gous to that of a micromaser. Less noisy than a laser, a 
micromaser [l^l can generate number-squeezed states 
of the cavity and exhibits not a single threshold transi- 
tion, but a series of transitions between different dynam- 
ical states. Although the SSET-resonator system and 
micromaser differ in the details of the interactions be- 
tween their respective sub-components, we find a num- 
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FIG. 1: (a) Schematic diagram of the SSET-resonator sys- 
tem. The SSET consists of a superconducting island linked 
by tunnel junctions to superconducting leads, across which a 
voltage Vds is applied, (b) Charge processes involved in the 
JQP current cycle (details are given in the text). 



ber of important similarities in their dynamics, many of 
which first arise when the resonator is sufficiently fast to 
match the time-scale of the electrical transport. Previ- 
ous theoretical studies of this system have concentrated 
on the limit of a slow resonator 0, 0] as it is only this 
limit which has so far been explored in experiments on 
nanomechanical resonators coupled to a SSET 2, J5] . The 
much faster resonator speeds which we also consider here 
might be achieved by making a smaller [ll| mechanical 
resonator. However, it should also be possible to use a 
superconducting resonator (e.g. a stripline resonator (l^ l 
for which much higher frequencies are practical. 

The SSET-resonator system we consider is shown 
schematically in Fig. [1^. A mechanical resonator acts as 
a voltage-gate with a position dependent capacitance and 
the coupling is controlled by varying the voltage applied 
to it (an analogous, though fixed, electrostatic coupling 
arises for a superconducting stripline resonator [l3|). The 
SSET is assumed to be biased close to a JQP reso- 
nance 0] where only three charge states of the SSET is- 
land are relevant. Current flows in a cycle (Fig.[T]D): co- 
herent Josephson tunneling between the left lead and the 
island produces a superposition of island charge-states |0) 
and 1 2), this is then followed by two quasiparticle tunnel- 
ing events between the island and the right lead which 
take the SSET island from state |2) to state |1) and fi- 
nally back to state |0). The resonator is modeled as a 
single-mode quantum harmonic oscillator with mass m 
and angular frequency lu. For a weak electrostatic in- 
teraction between the resonator and the SSET, the cou- 
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pling between the resonator position and the number of 
excess charges on the SSET island is Hnear 0, ■ The 
SSET-resonator system can therefore be described by the 
Hamiltonian 

E 

H = A£;|2)(2| - ^ (|0)(2| + |2)(0|) + ?iwa^a 

+Xsihu;'m/2y/^{a^ + a) (|1)(1| + 2|2)(2|) , (1) 

where a is the annihilation operator for the resonator, 
AE is the electrostatic energy difference between the is- 
land states |0) and |2), Xg is the displacement in the 
equilibrium position of the resonator when one electronic 
charge is added to the SSET island[§| and Ej the Joseph- 
son energy of the superconductors. The system can be 
considered to be in the weak-coupling limit when the volt- 
age applied across the SSET is much larger than the cou- 
pling energy, k = mto'^xl/eVds ^ 1 

The system evolves coherently under the action of the 
Hamiltonian, but dissipation arises from two sources: 
quasiparticle tunneling in the SSET and the resonator's 
surroundings. Thus the master equation for the SSET- 
resonator density matrix, p{t), derived using the Born- 
Markov approach 0, takes the form 

P = i^P= -J^[H,p\+ CleadsP + CdampmgP, (2) 

where the dissipative terms arising from the resonator's 
surroundings and from quasiparticle tunneling are de- 
scribed by the terms Cdamping and Cieads respectively. 
The resonator's surroundings are assumed to act like a 
thermal bath which wc describe using a Liouvillian (isj , 

CdampingP = " (f^ + 1) (a^ap + pfl^fl - 2apa)) 

which is guaranteed to preserve the positivity of the den- 
sity matrix. The temperature of the external bath is 
parameterised by the average number of quanta, n, that 
the resonator would have were it in thermal equilibrium 
with the bath and the resonator-bath coupling is given 
by a damping rate ^ext- The tunneling of quasiparticles 
from the island to the left lead is described by [8| 

AeadaP = - ^ [{ 1 2) (2 1 + 1 1) (1 1 , p} + 

-2(|l)(2| + |0)(l|)p(|2)(l| + |l)(0|)], 

where F is the quasiparticle tunneling rate. This is a 
simplified expression in which we have neglected differ- 
ences between the quasiparticle rates for the two pro- 
cesses, their variation with bias point and temperature 
dependence We have also neglected the de- 

pendence of the quasiparticle tunneling rates on the res- 
onator position as this is much less important than the 
coherent electro-mechanical coupling in the Hamiltonian 




FIG. 2: (Color online) Steady state properties of the resonator 
as a function of A_B, and resonator speed, uj/T: (a) average 
occupation number of the resonator, (n); (b) Fano factor of 
the resonator, both plotted on a natural-log scale, with k = 
0.01, 7er,t/r = 0.002, and n = Q. 

(Eq. [T|) d, [^, . These simplifications allow us to cap- 
ture the essential phenomenology of the system using a 
relatively compact model. Furthermore, simplified in this 
way our model of the SSET-resonator system close to the 
JQP resonance is dual to another system, that of a dou- 
ble quantum dot gated by a resonator ,15]. 

The value of /S.E determines the detuning of the 
Cooper pairs from resonance and can be changed con- 
tinuously by changing the applied gate voltage (this can 
be done independently of the coupling provided an ad- 
ditional, fixed, gate is usedQ). The sign of AE deter- 
mines the direction of (average) energy flow between the 
resonator and the SSET. For AE < the state |0) has 
more energy than |2), hence when a Cooper pair tun- 
nels onto the island it can pass some of its energy to the 
resonator (before quasiparticle decays occur), but when 
AE > the situation is reversed and hence the Cooper 
pair can absorb energy from the resonator. This means 
that when AE > the SSET damps the motion of the 
resonator, though because of the stochastic nature of the 
current, the resonator settles into a thermal- like steady 
state, not its ground state [!,[§]. However, when AE < 
the transfer of successive Cooper pairs passing from lead 
to island can pump the resonator and drive its oscillation. 

In order to explore the behavior of the system over a 
range of resonator speeds, from slow (w/F <C 1) to fast 
(w/F ^ 1) we solve the master equation using a numer- 
ical method. There are a number of methods which can 
be employed to calculate the steady state of a master 
equation. We have made use of the numerical routines 
implemented in the Quantum Optics Toolbox fl6| to cal- 
culate the steady state of the density matrix by evalu- 
ating the eigenfunction corresponding to the zero eigen- 
value of the Liouvillian in Eq. [2] (written in matrix form). 
Our method necessitates a truncation of the resonator's 
Hilbert space, which because of the need to ensure con- 
vergence, effectively translates into a lower limit on the 
value of the external damping we consider (we used a 
Fock state representation with up to 70 states). 

The effect on the steady state of the resonator of vary- 
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FIG. 3: (Color online) (a) Average occupation number and 
Fano factor of the resonator together with the most proba- 
ble number state in the density matrix, rimp, for uj/T — 0.1, 
Jext/T = 0.0005, K = 0.02 and n = 0. (b) Examples of the 
resonator Wigner functions: I fixed point state, II bistable 
and III limit cycle. The regions in (a) corresponding to the 
three basic Wigner function topologies are also marked as I 
(fixed point), II (bistable) and III (limit cycle). 



ing the detuning from resonance, AE, for a range of res- 
onators speeds, oj/T, is shown in Fig. 2 where the average 
number of resonator quanta, (n) (with n = a^a), and the 
Fano factor, F — {{n^) ~ are plotted. We as- 

sume (throughout) SSET parameters [9| F = Vds/eRj, 
Rj = h/e^ and Ej = hVds / {^&eRj) , where Rj is the 
SSET junction resistance. 

A more complete understanding of the resonator state 
is obtained from the Wigner transform of the reduced 
steady state density matrix of the resonator . Although 
the Wigner function of the resonator is different in de- 
tail for every set of parameters in Fig. 2, it turns out 
that because of the weak SSET-resonator coupling there 
are just three different topologies which arise: a state in 
which the resonator fluctuates about a fixed point, a limit 
cycle state in which the resonator undergoes finite ampli- 
tude oscillations and a bistable state in which both the 
fixed point and limit cycle type states coexist (illustrated 
in Fig. [3D. 

The behavior of the resonator is simplest in the limits 
of a very fast {(jj/T » 1) or very slow (w/F ^ 1) res- 
onator, where a wide separation of the SSET and res- 
onator time-scales limits their mutual interaction. In 
contrast, when lo/T ^ 1 the SSET and resonator inter- 
act most strongly and hence it is not surprising that it 
is in this regime that the most interesting features of the 
coupled dynamics first begin to emerge [17 1 . 

For uj/T > 1, Fig. 2 shows that the transfer of energy 
between SSET and resonator is increasingly concentrated 
around a series of lines. These lines are points where 
/S.E = ±jhuj with j an integer (the current through the 
SSET also shows peaks at these values). When AE' < 0, 
the resonator absorbs quanta of energy from the SSET 
leading to an enhancement of (n) and F. Because we 
are working in the weak coupling regime, the strongest 



effect occurs for AE = —Tiu, with progressively weaker 
signatures at AE = —2hu) and — 3?iw, as the exchange of 
multiple quanta result from higher order processes. A dip 
in (n) for AE — +fiuj is also visible though it becomes 
less pronounced as oj/T is increased. This dip occurs, 
despite the fact that we have assumed an external bath 
at zero temperature, as the noise arising from coupling 
to the SSET excites the resonator out of its ground state 
even for AE > 0. These results are consistent with cal- 
culations for a resonator coupled to a double quantum 
dot in this regime [l^l ■ 

In the limit of a slow resonator, w/F <C 1, variation of 
AE leads to two distinct transitions in the state of the 
resonator as is shown in Fig. [3] The transitions become 
sharper as the ratio of the SSET current to the relax- 
ation rate jext increases, which suggests that they can 
be thought of as non-equilibrium phase transitions in a 
'thermodynamic limit' where this ratio diverges [lol. ITst. 
The most probable number state in the density matrix, 
n-mp, provides a convenient order parameter of the sys- 
tem. For AE/eVds ^ 0.4, the resonator is very close 
to being in a thermal state and hence there is a single 
peak in the Wigner function, but as AE is decreased 
and approaches zero the resonator state broadens (devi- 
ating strongly from a thermal state), until for AE ~ 
a limit cycle begins to appear. The transition to a limit 
cycle is continuous as the radius of the cycle (and n,„p) 
grows steadily from zero as AE is decreased. However, 
eventually (as AE becomes suSiciently negative) the sys- 
tem passes through a region of bistability where both the 
limit cycle and a central peak co-exist in the Wigner func- 
tion before the limit cycle disappears entirely and only 
the central peak remains. This behavior implies a dis- 
continuous transition and is marked by a jump in rimp- 

For a very slow resonator, the current flowing through 
the SSET decays monotonically to zero far from the cen- 
ter of the JQP resonance. This explains why the limit 
cycle eventually disappears as AE is decreased: the en- 
ergy per unit time transferred to the resonator is propor- 
tional to the current and as this decays, the system will 
eventually be stabilized by the external damping. 

A striking feature of the intermediate regime (w/F ~ 
1) in Fig. 2 is the drop in the Fano factor of the res- 
onator which occurs in the limit cycle state. Remark- 
ably, there is a region where the Fano factor falls below 
unity, implying that the resonator is driven into a number 
squeezed (i.e. non-classical) steady state. We attribute 
this squeezing to non-linear damping induced by the in- 
teraction with the SSET . Squeezing of this kind is a 
characteristic of the micromaser ^, ilOi] , but never occurs 
in the (more noisy) conventional laser [l3| . As with the 
micromaser, we find that increasing the noise in the sys- 
tem (by increasing the temperature of the bath) rapidly 
washes out the squeezing eff'ect. 

The regime where w/F ~ 1 is also where further tran- 
sitions beyond a simple limit cycle state first occur as 
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FIG. 4: (Color online) Changes in the resonator state (a) 
and corresponding Wigner functions (b) as the coupling k is 
increased for AE/eVds = -0.1, tu/F = 1, 7e^t/r = 0.001 and 
n = 0. After an initial continuous transition {k ~ 0.025^), the 
system undergoes discontinuous transitions at k ~ 0.22^^ and 
K ~ 0.44^. The inset in (a) shows the region where the Fano 
factor drops below unity. From top to bottom, the Wigner 
functions show: a fixed point state, a single limit cycle and 
two met ast able limit cycles. 

the SSET-resonator coupling is made stronger. A hall- 
mark of the micromaser is the sequence of transitions 
which occur as the atom-cavity coupling is increased [lo| . 
We find that the resonator undergoes a very similar se- 
quence of transitions as the coupling is increased at fixed 
AE{< 0), the first of which is continuous (for \AE\ suf- 
ficiently small) and the rest discontinuous. An example 
of this behavior is shown in Fig. H) As the coupling is 
increased from zero the resonator first undergoes a con- 
tinuous transition between fixed point and limit cycle 
states. Within the limit cycle state F drops below unity, 
but as K increases above 0.01, F grows sharply as a sec- 
ond metastable limit-cycle state begins to appear. For 
K ^ 0.04, the Wigner functions become more compli- 
cated, showing further metastable limit-cycles together 
with regions where parts of the Wigner function have neg- 
ative values (indicating non-classicality of the resonator 
state). The relative weights of the different limit cycle 
states change with k, leading to discontinuous transitions 
signalled by jumps in Ump as the most probable state of 
the resonator changes. In Fig. 2] only the first (continu- 
ous) transition shows sharp features in (n) and F as we 
are far from the 'thermodynamic limit' [lO| (we have cho- 
sen a relatively large jext to ensure that our truncation 
of the oscillator state space remains valid). 

In conclusion, we have analyzed the dynamics of a res- 
onator coupled to a SSET near the JQP resonance. The 
SSET-resonator system has a rich dynamics with many 
features in common with the micromaser including con- 
tinuous and discontinuous transitions between resonator 
states and intrinsically quantum features such as non- 



classical steady-states. In practice, the resonator dynam- 
ics could be inferred from signatures in the current and 
current noise [l^ of the SSET or, for a superconducting 
resonator, by coupling it to a transmission line [l^ . The 
effects described are strongest when uj/T ^ 1 and will 
require minimal thermal excitation (n < 1) for experi- 
mental observation of the non-classical features, implying 
high resonator frequencies. Similarities with the micro- 
maser suggest the possibility of observing critical slowing 
down, hysteresis and quantum jumps in the resonator dy- 
namics [Gj . 

After submission of this paper, [2^ was published 
which explores analogies between nano-electromechanical 
systems and laser physics in the slow-resonator limit. We 
thank T. Brandes, M. Blencowe and T. Harvey for useful 
discussions and acknowledge funding from the EPSRC. 



[1] M.P. Blencowe, Phys. Rep. 395, 159 (2004). 

[2] R.S. Knobel and A.N. Cleland, Nature 424, 291 (2003); 

M.D. LaHaye et al.. Science 304, 74 (2004). 
[3] A. Naik et al.. Nature 443, 193 (2006). 
[4] A.A. Clerk and S. Bennett, New J. Phys. 7, 238 (2005). 
[5] T. Novotny, A. Donarini and A.-P. Jauho, Phys. Rev. 

Lett. 90, 256801 (2003); D. Fedorets et al, Phys. Rev. 

Lett. 92 166801 (2004); D.W. Utami et al., Phys. Rev. B 

74 014303 (2006); O. Usmani, Ya.M. Blanter and Y.V. 

Nazarov, /cond-mat/0603017 (unpublished). 
[6] H. Wahher et al.. Rep. Prog. Phys. 69, 1325 (2006). 
[7] I. Bargatin and M. Roukes, Phys. Rev. Lett. 91, 138302 

(2003). 

[8] M.S. Choi et al., Europhys. Lett. 53, 251 (2001). 
[9] M.P. Blencowe, J. Imbers and A.D. Armour, New J. 
Phys. 7, 236 (2005). 

[10] P. Filipowicz, J. Javanainen and P. Meystre, Phys. Rev. 
A 34, 3077 (1986). 

[11] X.M.H. Huang et al.. Nature 421, 496 (2003). 

[12] A. Blais et al., Phys. Rev. A 69 062320 (2004); A. Wall- 
raff et al.. Nature 431, 162 (2004). 

[13] C. Flindt, T. Novotny and A.-P. Jauho, Phys. Rev. B 70, 
205334 (2004). 

[14] The numerical method can also be used to explore the be- 
havior of a more accurate (though less general) model of 
the SSET resonator system. Indeed, we have checked that 
including the position dependence of the quasiparticle 
rates does not qualitatively alter the results we present. 

[15] T. Brandes and N. Lambert, Phys. Rev. B 67, 125323 
(2003). 

[16] S. M. Tan, J. Opt. B: Quantum Semiclassical Opt. 1, 424 
(1999). 

[17] D.A. Rodrigues and A.D. Armour, New J. Phys. 7, 251 

(2005) . 

[18] H. Haken, Rev. Mod. Phys. 47, 67 (1975). 
[19] C. Flindt, T. Novotny, and A.-P. Jauho, Europhys. Lett. 
69, 475 (2005). 

[20] S. Bennett and A.A. Clerk, Phys. Rev. B 74, 201301 

(2006) . 



